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It is widely believed that integer quantum Hall systems do not have fractional excitations. Here 
we show the converse to be true for a class of systems where integer quantum Hall effect emerges 
spontaneously due to the interplay of itinerant electrons and noncoplanar magnetic ordering. We 
show that magnetic Z2 vortices in such systems, whose topological stability is guaranteed by the 
structure of the order parameter, indeed carry fractional charge. 



The discoveries of the integer and fractional Quantum 
Hall (QH) effects have revolutionized condensed matter 
physics: the important concept of a topological invariant 
was introduced to explain the quantized Hall conductiv- 
ity of the former [1], while the novel notion of topologi- 
cal order, i.e., a type of nonlocal order with no Landau 
symmetry-breaking and no local order parameter, was in- 
troduced to describe the latter [2] . Topological order goes 
hand in hand with exotic phenomena such as fractional 
charge and statistics [2]. While the strongly-correlated, 
topologically-ordered fractional QH systems indeed have 
fractional quasi-particles [3, 4], their more traditional in- 
teger counterparts are assumed to only support natural 
excitations with integer charge. Contrary to this assump- 
tion, we predict in this paper that factional excitations 
can generically emerge in a class of integer QH systems. 

The existence of fractional excitations in a fractional 
QH liquid follows from the Laughlin's argument: upon 
adiabatic local "insertion" of a flux quantum, a fractional 
charge q = a xy e, where a xy is the fractional Hall conduc- 
tivity, flows in from infinity. Since integer flux quantum 
can be "gauged away," the charge q is in fact the charge 
of an elementary quasiparticle. Naturally, insertion of a 
fractional flux can also lead to the same result with an 
integer Hall response. Fractional fluxes, however, cannot 
be gauged away, and do not naturally appear in a tradi- 
tional integer QH system, i.e., two-dimensional electron 
gas in a magnetic field. (In principle, half-vortices could 
be introduced in artificial heterostructures, by placing 
a QH system in proximity to a superconductor [5].) It 
may thus appear that the presence of intrinsic fractional 
excitations in a QH system requires a fractional a xy . 

However, besides topological order, charge fractional- 
ization can also occur through another mechanism: cou- 
pling of noninteracting electrons to topologically stable 
defects in ordered backgrounds, such as a pattern of lat- 
tice distortions in polyacetylene or graphene-like struc- 
tures [6-9] or in a superconducting vortex [10-13]. The 
presence of an ordered background requires spontaneous 



symmetry breaking, which does not happen in traditional 
integer QH systems. But QH effect can also emerge as 
a result of spontaneous time-reversal symmetry breaking 
in the absence of a magnetic field [14-18, 20]. As we will 
show in this paper, such spontaneous integer QH systems 
can support topologically stable defects, which remark- 
ably act as an effective fractional flux, thus giving rise to 
natural excitations with fractional charge. 




FIG. 1: (color online), (a) Chiral spin ordering on a triangular 
lattice. Four orientations of the local magnetic moments cor- 
respond to the normals to the sides of a regular tetrahedron, 
(b) Vortex configuration. 

Spontaneous QH effect can be induced by noncopla- 
nar magnetic ordering, as depicted in Fig. la, serving 
as the background for itinerant electrons [16-18]. Even 
though in this work we focus on the specific model of 
Ref. [18] on the triangular lattice, we would like to em- 
phasize that our qualitative conclusions are independent 
of the microscopic details, and rely on only two ingredi- 
ents: noncoplanar magnetic ordering and the presence of 
an incompressible integer QH electronic state. We can, 
for instance, apply the same arguments to the model of 
Ref. [16] on the kagome lattice, as shown in the supple- 
mentary material [19]. 

In the model of Ref. [18] (Fig. la), the magnetic order- 
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ing spontaneously and completely breaks 0(3) symmetry. 
The existence of topologically stable vortex defects in a 
noncoplanar magnetically ordered medium follows from 
the nontrivial fundamental homotopy group of the order 
parameter space, tti(0(3)) = m(SO(3)) = Z 2 [19, 24]. 
An example of the vortex spin texture is shown in Fig. 
lb. For instance, it can be obtained by rotating the or- 
der parameter (every magnetic moment) by a position- 
dependent angle 

<j)(r) = vaxg(x + iy) (1) 

around the z axis (assuming the vortex core is at the ori- 
gin) ; v is the winding of the vortex. The two topological 
classes are set by the parity of v. Unless stated otherwise, 
here we consider the v = 1 case. 

A model that gives rise to the spontaneous QH state 
is the following Kondo lattice model: 

H = -UjC^Cja - Hc\ a Cia + Jc|q S » • ValiCiji, (2) 

where electrons hop on a two-dimensional (2D) lattice 
and interact with the local magnetic moments via on- 
site exchange interaction. Summation over repeated site 
(roman) and spin (greek) indices is implied. Here tij is 
the intersite hopping, J is the exchange interaction con- 
stant, Si is the local magnetic moment (assumed classi- 
cal), <r = (a x , a y , <t z ) is the vector of Pauli matrices, and 
c ia is the operator of electron annihilation on site i with 
spin a. On a triangular lattice, the model has a spon- 
taneous QH groundstate (Fig. la.) for several electron 
densities [18, 21-23]. 

The magnetic texture may become inhomogeneous ei- 
ther due to thermal or quantum fluctuations, or due 
to the presence of topologically stable defects, such as 
vortices (Fig. lb). Such inhomogeneous states can be 
mapped onto a state with homogeneous order parame- 
ter, but in the presence of an effective (in general) non- 
Abelian gauge potential. This mapping allows one to 
conveniently calculate the charge and spin currents in re- 
sponse to the order parameter distortions. In particular, 
the vortex configuration corresponds to a spatially local- 
ized non-Abelian flux, which can be used to determine 
the charge of the vortex. Indeed, an inhomogeneity cor- 
responds to a rigid position- and possibly time-dependent 
distortion of some reference state, Si = 1Z(ri,t)S®. The 
rotation of the order parameter in the classical spin space 
can be transformed into a unitary rotation U(r i} t) = Ui 
of the electron spinors, according to U\a ■ SiUi = er ■ S°. 
Introducing new fermions ipi — U(ri,tyCi (spin index is 
suppressed), the Hamiltonian (2) becomes 

h= -i^jujdtU^-u^juju^ 



+ JS° ■ 4<"l>i. (3) 

Assuming that the variation of the texture is slow on 
the lattice constant scale, we can make an expansion, 
Ujllj = Uj[Ui + {vj - n) ■ Vt/j]. It is convenient to in- 
troduce SU(2) vector potential A" = -iWd"U = A v a a a , 
with the indices v — {t, x, y} representing the space-time 
components and a = {1, 2, 3} the SU(2) generators. The 
Hamiltonian can then be written as 

H = H - r a A» a , (4) 

where Hq is the Hamiltonian corresponding to the static 
uniform spin structure, and the currents are defined as 

J a ° - V(ri,t)W(ri,t), (5) 

J a = iUj{Xi-Xj)l()(Ti,t)^(Ta^{rj,t), (6) 

J a = i-t l j{y l ~y j )^{-r l ,t)^a a tp{Y jl t). (7) 

Denoting the 2x2 unit matrix by ao, for a = the 
definitions above also give the physical charge density 
and current operators [32]. 

If the texture is slowly varying in time (compared with 
the inverse energy gap in the spectrum), as well as in 
space, the vector potential A v a is small and the current 
expectation values can be calculated with the linear- 
response theory. Through an explicit calculation for the 
model of Ref. [18], we find the following nonvanishing 
responses [19]: 

a° x ° y = ~e 2 /h, a a x ^e 2 /3h (for a ^ 0), (8) 

which determine the charge and spin currents through 
(J%) = a% a y e r)liv d 11 A v a . The signs of both conductivities 
are flipped by switching between 3 /4 and 1 / \ fillings, or 
by changing the sign of the chiral ordering. 

The vortex texture is necessarily singular near the core; 
therefore, the gauge transformation that one needs to ap- 
ply to unwind it will be singular as well. The simplest 
transformation that takes the vortex texture into a uni- 
form one is e la 3<t>( r )/' 2 j where <j> is the angle of rotation (1) 
around the z axis. However, since upon going around the 
vortex, 4>(y) — > <j>(r) + 2tt, the unitary changes sign, this 
would correspond to antiperiodic boundary conditions for 
fermions along a line connecting the vortex to infinity. To 
avoid this complication, the above SU (2) transformation 
can be augmented by a £7(1) one [25]. The combined 
transformation U(r) = e*^ <73+1 ^^ r ^ 2 is only acting on 
up-fermions. Its associated gauge potential is 



which has the field strength zero everywhere except for 
the vortex core. 
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Due to the singular nature of the gauge potential, we 
cannot directly apply the linear-response formalism in 
the vicinity of the vortex core. To calculate the vortex 
quantum numbers, we can instead invoke an analogue 
of the Laughlin argument. The flux of the non-Abelian 
gauge field through the vortex core is 

* = <j> Adv = (1 + 03)71". 

Now, suppose that the flux is turned on adiabatically 
from zero to <I>. That will generate a non-Abelian cmf 
acting on electrons, which at large distances from the core 
will be nearly uniform (tangential to any circle centered 
at the vortex). The vortex quantum numbers are then 
obtained by integrating the associated currents generated 
in response to this emf. 

Before proceeding, let us note that the Laughlin argu- 
ment relies on two conditions. First, we need to have a 
continuous sequence of gapped Hamiltonians connecting 
the one with flux zero to the one with <&. Second, we 
need a continuity equation relating the currents we can 
calculate in linear response to quantum number densi- 
ties. We have explicitly identified a sequence of gapped 
Hamiltonians in the limit of large J and we thus expect 
that an adiabatic process exists for arbitrary J. The 
quantum numbers of interest for us are charge and spin. 
Since total electron number commutes with the Hamil- 
tonian, the charge current strictly satisfies a continuity 
equation. We do not have such continuity equation for 
the spin current, and thus our vortex defects do not have 
a well-defined spin quantum number. 

However, since the divergence of the induced current 
is zero far from the vortex, it is expected that the spin 
current will be nearly conserved, except near the vortex 
where it accumulates. This suggests that the magnetiza- 
tion attached to a vortex might be still close to the ex- 
pected value obtained from integrating the spin current. 
With the assumptions above, since the final flux both in 
spin-(73 and charge channels is half of the flux quantum, 
the accumulated charge and magnetization should be [33] 

q = a° x ° y /2, m z « afjl. (9) 

(the extra 1 / 2 in the expression for m z is due to the fact 
that electron spin is tr/2.) For less symmetric systems 
such as in Ref. [16] there are contributions from mixed 
conductivities, like a® y or for instance, to q and m z 
as well as to other spin components. 

There are also other possible choices of the gauge trans- 
formation that unwinds the vortex, e.g., with the gauge 
potential A v = (-1 + a 3 ) d u <f)/2. While such different 



choice does not change the magnetization, the charge ac- 
cumulated in the vortex core changes sign. This am- 
biguity is naturally understood in terms of the electron 
occupancy of a particular localized electronic state, e , 
inside the spectral gap. When this state is empty the 
charge of the system is — 1 / 2 and when it is occupied it is 
+ 1 /2, all relative to the uniform state. In general, there 
can be more than one localized state inside the vortex 
core. Occupying any of these states increases the vortex 
charge by one electron charge [this corresponds to more 
general gauge choices, K v = (I + 2n + 03) d v <fi/2, with n 
any integer]. However, this cannot change the fact that 
the charge of the vortex has to remain half-odd-intcgcr 
for an odd vorticity. On the other hand, for an even 
vorticity, the charge induced according to the Laughlin 
argument will be integer. This is consistent with the ho- 
motopy classification that says that double vortex can be 
smoothly connected to a uniform state. 

We now check the above results numerically using the 
explicit model of Ref. [18] (Fig. I). We plot, in Fig. 2, 
the charge and magnetization distribution in the vicin- 
ity of the vortex core. As expected, the charge local- 
ized in the core is half-odd-integer for odd winding and 
integer for even winding. The agreement between the 
vortex magnetization obtained numerically with what is 
expected from Eq. (9) is not perfect because the spin cur- 
rent is not conserved. Nevertheless, particularly for large 
J, the discrepancy is not too large and we verify that the 
vortex spin polarization approximately scales with the 
vortex winding, as shown in Fig. 3b. 

We have also considered deviations from the fully sym- 
metric assumptions. In particular, we added a Zee- 
man field h along z axis acting on electrons, i.e H — > 
H + hJ2i c ia a a0 c ii3- We verified that as long as the spec- 
tral gap does not close, the charge Hall conductivity does 
not change. However, in contrast to the fully symmet- 
ric case, new nonzero response functions emerge, namely 
a° y , a y x 7^ 0, which correspond to the charge response to 
A3, the CT3 component of the gauge potential. In Fig. 3a, 
we plot the dependence of the vortex charge on h. The 
solid line is q = —a xy /2 — (a®y — (Ty X )/4 — 1, which di- 
rectly follows from the application of the Laughlin argu- 
ment. Note that the charge is generally irrational and 
determined modulo an integer. The agreement is very 
good, all the way to the value of h where the gap in the 
electronic spectrum closes. 

The energetics of Z 2 vortices is similar to the usual Z 
vortices for SO(2) order parameter [24]: i.e., the energy 
of an isolated vortex scales logarithmically with the sys- 
tem size. Nevertheless, pairs of log-confined vortex pairs 
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FIG. 2: (color online). Charge (a) and magnetization (b) 
density distributions around a vortex. The vortex is indicated 
by the green dot at the center. The system parameters are as 
follows: J — > oo, L — 30, h = at 3 / \ electronic filling. 



(a) J/t = 4.0, L = 40, v = 1 (b) J/t — > oo, L = 30, ft = 
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FIG. 3: (color online), (a) Zeeman field dependence of the 
vortex charge. Blue circles were obtained by exact diago- 
nalization. The red solid line is the expected result from a 
Laughlin adiabatic argument using the uniform dc Hall con- 
ductivities of an infinite uniform system shown in Fig. la. 
(b) Net magnetization accumulated nearby the vortex for dif- 
ferent windings and with rotations around the 3 axes. The 
solid line corresponds to Eq. (9). Both results were obtained 
for / 4, electronic filling. 

will appear due to thermal fluctuations at finite temper- 
atures. In addition, inclusion of quantum spin dynamics 
may lead to an intriguing possibility that the quantum 
fluctuations transform the noncoplanar ordered state into 
a chiral spin liquid [26-28] at zero temperature. There, 
the fractionally charged Z 2 vortices discussed above may 
turn into deconfined point-like excitations. 

Finally, we comment on the vortex exchange statistics. 
There are two distinct possibilities for the case of two 
vortices: (1) the total charge of two vortices is even, and 
(2) the total charge is odd. The former case will be re- 
alized if the vortices are pulled apart from the uniform 
"vacuum": since the initial state has total charge zero 
(relative to the uniform background), the state with two 
vortices will keep the same charge. Since the charge of 
an individual vortex is half-odd-integer, for large inter- 
vortex separation there are two energetically equivalent 
ground states that correspond to vortex charge config- 
urations (V2, — V2) an d (— V27 1 /2)- This degeneracy 
can lead to non-Abelian exchange statistics, but unlike 
the Majorana states in the p-wave superconductor-like 
systems [12, 13], there is no topological protection. In 
other words, any local disorder can shift the bound state 



energy in a given vortex and lift the degeneracy. 

The case (2) can be obtained, e.g., upon electron or 
hole topological doping of the system with equally (half- 
odd-integer) charged vortices. In this case, there is no 
ground-state degeneracy and exchanging two vortices can 
only give rise to an Abelian Berry phase. The anyonic 
character of the vortices follows from a simple argument. 
Consider two indistinguishable distant vortices with sin- 
gle winding and equal charge. Near each vortex the gauge 
potential is A" — (1 + 03) d"4>/2, with <j> being the angu- 
lar coordinate relative to the given vortex core. There- 
fore, each vortex has a non-Abelian gauge flux attached 
to it, equivalent to a U(l) flux 2-7T in the spin-up chan- 
nel. The number of up electrons occupying each vortex 
is n-f = (q + 2m 2 )/2. The statistical angle due to the 
flux attachment is the product of the charge and half 
of the flux, which for the idealized case of Eq. (9) is 
9 = tt(<7°° +(T^)/4. Note that the deviations of m z from 
this idealized case are small for J — > 00 (see v = 1 data 
in Fig. 3b) . The numerical observation of the statistical 
angle, however, is complicated by the spin nonconserva- 
tion. 

In summary, we have shown that spontaneous integer 
quantum Hall systems, which emerge from the interplay 
of itinerant electrons and noncoplanar magnetic ordering, 
generically support topologically stable excitations with 
fractional charge and anyonic statistics. Promising can- 
didates materials which may exhibit this physics could 
be the systems of Na^CoC^ type, which near x = 0.5 are 
known to have a noncollincar order, as well as anoma- 
lously large Hall response [31]. The fractional charge 
predicted in this work may be accessible through direct 
imaging of the local charge profile, as shown in Fig. 2, 
e.g., by scanning force microscopy. Also anyonic ex- 
change statistics may have unusual consequences in real 
materials. Perhaps the most intriguing among them is 
the possibility of the Anyonic Superconductivity [28-30] . 
When a system is doped away from the chiral Mott in- 
sulating state, it may energetically prefer to accommo- 
date the carriers by creating vortices with intragap states. 
As we have just argued, such occupied vortex states are 
anyons, which, at finite density and low enough temper- 
ature, may go into a superconducting state [30]. 
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I. 50(3) VORTICES 

A rotation in 3D can be parameterized by an angle and an axis. Thus, the order-parameter space 
can be geometrically represented by a solid sphere of radius tt with antipodal points on the surface 
identified: the distance from each point to the center of the sphere represents the angle of rotation, 
while the vector connecting the point to the center gives the axis of rotation. The identification of 
antipodal points follows from the fact that clockwise and counterclockwise rotations by angle tt are 
equivalent. A ID loop in real space maps onto a ID loop in the order-parameter space. As seen in 
Fig. la, due to th identification of antipodal points, there are two topologically distinct loops in 
the order-parameter space: contractible (topologically trivial), and noncontractible (topologically 
notrivial). The noncontractible loop corresponds to a nontrivial vortex. For example, the vortex 
configuration described in the main text corresponds to the following noncontractible loop: a 
straight line passing through the North pole, the center of the sphere and the South pole (which 
is identified with the North pole) of Fig. lb. 




FIG. 1: a) The order parameter space and the two types of topologically distinct loops: a trivial loop and a 
vortex. The loops can lie on the surface of (black) or inside (gray) the solid sphere. The loops corresponding 
to a vortex connect two antipodal points on the surface of the sphere, b) The nontrivial loop corresponding 
to the vortex shown in the main text. 



II. CALCULATION OF THE RESPONSE FUNCTIONS 

In this section, we present some details regarding the calculation of different response functions. 
We use the general linear-response relation for the expectation value of different currents: 

(r a (r,t))= [ R2(t-t',v-v')Al(t',v% (1) 

Jt',r' 

where FfV are the linear-response functions calculated with respect to the uniform state. The 
number of non-zero response functions can be quite small due to symmetries. In any QH system, 
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the relationship between the charge current and the electromagnetic potential Aq is given by 
(Jq) = cr^e^i,^ Aq, where a®® is just the regular charge Hall conductivity. Comparing with 
Eq. (1), this gives 

lQ{t -t',r- r') = af y 5{t - t')5(r - v')e v ^. 

Other response functions R^it — t') = —i([J^(t), J^(t')])9(t — t') involve either spin-spin or spin- 
charge correlators. For the particular model of Ref. 1, the only other nonvanishing Hall response 
is in the spin-spin sector, namely a^y, a = 1 . . .3. That is due to the symmetric nature of the 
ordered state: a combination of a spatial translation and a spin rotation is a symmetry. For other 
less symmetrical systems, such as an arbitrary spin configuration in the model of Ref. 2, the Hall 
conductivities will have a less structured form. 

All conductivities can be computed from the general expression 

2 i I ix\nmk( rJAmrtk 

<=h2iTi £ - Me mk )]. (2) 

n,m,k 

The expression above is applicable to translationally invariant systems where momentum k is a 
good quantum number. For a unit cell of M sites, the momentum-space Hamiltonian %(k) can 
be generically written as a 2M x 2M matrix. (The factor of two accounts for electron spin.) For 
each momentum k, we can diagonalize this 2M x 2M Hamiltonian and obtain its eigenvalues and 
eigenvectors: %(k) = e m k|?rt(k)), m = 1 . . . 2M. The eigenvalues e m k give the 2M energy bands. If 
we have symmetries, as in the triangular lattice case discussed below, these bands can be degenerate. 
The matrix elements (J) nmk = (n(k)| J(k)|m(k)), which appear in Eq. (2), can be constructed 
explicitly using the eigenvectors \m(k)) and |n(k)), where J(k) is an appropriate current operator 
written as a 2M x 2M matrix in the same basis as "H(k). Also, np(e) = (l + e 13 ^^) is the 
Fermi occupation number at a given chemical potential [i and inverse temperature /3. Throughout 
this work, we focus on zero temperature, where np is a step function. 

The only ingredients for computing Eq. (2) are then the 2M x 2M Hamiltonian H(k) and 
the corresponding 2M x 2M charge and spin current operators J% ,y (k) for a = . . . 3. With 
these ingredients, one can diagonalize "H(k) to obtain the eigenvalues and eigenvectors, use the 
eigenvectors to construct the matrix elements of the current operators, perform the sum over m 
and n, and finally integrate the resulting expression over momenta k in the Brillouin zone. 



III. TRIANGULAR LATTICE 

In this section, we present the details of the calculation of the response functions [Eq. (9) of 
the main text] on the triangular lattice. Based on the above discussion, our main task is to write 
"H(k) and Ja ,y (k)- We first choose an explicit tetrahedral magnetic ordering represented by the 
moments S a , a = 1 . . . 4 in Fig. 2. With four sublattices and two spin species, we can write the 
Hamiltonian as an 8 x 8 matrix in momentum space. As in the main text, we assume an additional 
Zeeman field h in the z direction. We choose the following basis: 

= ( C lfk' C Uk' 4tk' C 24.k' C 3fk' C 34-k' C 4fk' C li-k)' 
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to write the Hamiltonian, i.e., 



^T = ^* t k ^ T (k)* kj 



where the subscript T indicates the triangular lattice. We can then write 



ttr(k) = J 



( Si ■ a + ha 3 

S 2 ■ a- + ha 3 



\ 



where 



£r(k) 








/ o 

— 2t cos(k • ai) 
— 2t cos(k • as) 
y — 2t cos(k • a2) 








\ 



S3 • cr + /io- 3 

S^-cr + ha 3 J 



+ £ T (k) 8) a c 



-2tcos(k-ai) — 2tcos(k-a3) — 2t cos(k • a2) \ 

— 2tcos(k-a2) — 2tcos(k-a3) 

-2icos(k-a2) — 2icos(k-ai) 

-2tcos(k-a 3 ) -2tcos(k-ai) / 
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FIG. 2: An explicit chiral configuration of magnetic moments with tetrahedral ordering on the triangular 
lattice. For site i in sublattice a = 1 . . . 4, S» = S a as shown in the figure. The vectors are the lattice 
vectors. 

Charge (a = 0) and spin (a = 1,2,3) current operators can then be simply written in the same 
basis as 



J*(k) = d k J T (k) a\ jy(k) = d k £ T (k) a a . 



(5) 



IV. K AGO ME LATTICE 



As stated in the main text, the mechanism for fractionalization is generic, and does not depend 
on the specifics of the triangular lattice model. To demonstrate this fact, we consider another 
system where noncoplanar magnetic ordering leads to spontaneous integer quantum Hall effect. 2 
The model is on the kagome lattice, i.e., a lattice of corner-sharing triangles, shown in Fig. 3. 
In Ref. [2], a limit of large Kondo coupling was assumed, where the problem reduces to one of 
spinless electrons hopping on the lattice with a magnetic flux <f> (—2(f)) inserted in each triangular 
(hexagonal) plaquette. We choose the following noncoplanar ordered texture represented in Fig. 3: 
the projections of the magnetic moments S a , a = 1 ... 3, on the xy plane form 120-degree ordering 
on this plane. All the moments, however, are canted down perpendicular to the xy plane, creating 
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a noncoplanar texture. In the limit of large J, this setup reduces to the model of Ref. [2]. For con- 
creteness, we work with the specific moment configurartion shown in Fig. 3, where three moments 
point to the corners of a tetrahedron, as in the case of the triangular lattice. Other canting angles 
would give rise to similar physics. In analogy with the triangular lattice case in the main text, we 
introduce as an extra tuning parameter in the Hamiltonian a Zeeman field h in the z direction. 

First, we calculate the different response functions a®y and a^, for a = ... 3, which can lead to 
the accumulation of charge, when adiabatically turning on a Z 2 vortex, via the Laughlin argument. 
Then, we compute the charge accumulation around the vortex by exact diagonalization in a finite 
system. We obtain very good agreement between the two methods. Interestingly, in the absence 
of a Zeeman field, we obtain a charge 1/2, just as in the traiangular lattice case, when the spin 
rotation axis lies in the xy plane. 

As before, to calculate the responses, we simply need to write down %(k) and Ja ,y (k). 




51 = (-V5/V3, -V2/3, -1/3) 

5 2 = (V2/A-V2/3.-1/3) 

5 3 = (0,2^/3,-1/3) 



FIG. 3: An explicit chiral configuration of magnetic moments on the kagome lattice. For site i in sublatticc 
a = 1 . . . 3, Sj = S a as shown in the figure. The vectors are the lattice vectors. 



Defining 



t J J J J 



V c ltk' C 14k' C 2tk> u 24.k> °3tk' u 34.k 



c 3tk J c : 



we can write the Hamiltonian as Hk = Y2k ^k^fOO^kj where the 6x6 'HR-(k) matrix is given 
by 



/ 5i • <t + ha 3 



U K (k) = J 



V 





S 2 -(r + ha 3 I +£ K (k)®a°, 

S 3 ■ a + ha 3 



(6) 



with cr° representing the 2x2 identity matrix and the 3x3 matrix £ft:(k) given by 



( ° 

— 2t cos(k • ai) 

l — 2t cos(k • a2) 



-2t cos(k • ai) 


-2t cos(k • 83) 



-2t cos(k • a2) 
-2t cos(k • as) 




(7) 



The lattice vectors a^ are defined in Fig. 3, and the subscript K indicates the kagome lattice. The 
spectrum of the Hamiltonian above consists of 6 bands, as shown in Fig. 4, e.g., for J = 6 and 
h = 0. In the limit of large J, the three lower bands are separated from the three upper ones by 
an energy scale proportional to J. Throughout this section, we only consider the 5/6 filling, i.e., 
we work at a chemical potential where only the highest energy band is unoccupied. This is one 
of the several fillings where the system exhibits integer quantum Hall response 2 , as seen in Fig. 4. 
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FIG. 4: The band structure of Hamiltonian (6) for J = 6 and h = 0, where k = fciQi + &2Q2, and 
Qi = (1, — l/y3) and Q2 = (0,2/v3) are reciprocal vectors. The Chern numbers of the bands are shown 
on the right-hand side. From the Chern numbers, we find that the system exhibits integer quantum Hall 
effect at 1/6, 2/6, 4/6 and 5/6 filling fraction. 



Charge (a = 0) and spin (a = 1,2,3) current operators can then be simply written as in Eq. (5), 
with St replaced by Sk- The response functions for h = have the following properties: 



T oo _ 00 _ 2 /h 01 ra 01 q 02 02 03 _ 03 

xy v x ' ' y£ 5 yxi xy yx' 



As expected all the response functions saturate to fixed values for large J, where the model reduces 
to the model of Ref. [2]. 
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FIG. 5: a) Vortex charge for the kagome lattice with different coupling strengths. Red (blue) dots were 
obtained by exact diagonalization for a vortex with rotations around the z (x) axis. The solid lines are the 
analytical predictions from Lauglin's argument argument. Data for y-axis vortices are very similar to the x 
one and are not shown for clarity, b) Zeeman-field dependence of the vortex charge for the kagome lattice. 
Red dots were obtained by exact diagonalization. The solid line is the result expected from Laughlin's 
adiabatic argument. The electronic filling is 5/6 in both cases. 



Following the argument in the main text, the charge trapped by a vortex with the spin rotation 
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axis around the direction n = x, y, z is then given by 

^ = ^["</2-«-^)/4]. (8) 

Fig. 5a shows the charge of a vortex for several coupling strengths. The horizontal axis does 
not have a linear scale so both the variations at small J and the saturation for large J can be 
displayed. The results computed from Laughlin's argument show very good agreement with the 
numerical ones. Fig. 5b shows the charge q~ z trapped by a vortex as a function of the Zeeman field 
h along z axis. Like the triangular lattice case, we find that the charge Hall conductivity cr^ does 
not change as long as the spectral gap remain open. The off-diagonal responses a®? y on the other 
hand change as a function of h resulting in a continuously changing, generally irrational, fractional 
charge. The results computed from the Laughlin's argument again show excellent agreement with 
the ones obtained by numerical diagonalization on finite lattices of size 24 x 24. 
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